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I.  INTRODUCTION 

'~^The  spectroscopic  study  of  van  der  Waals  (vdW)  clusters  in  a 
molecular  jet  supersonic  expansion  has  greatly  increased  our 
understanding  of  these  species.  Additionally,  clusters  composed  of  an 
aromatic  molecule  (solute)  surrounded  by  small  hydrocarbons  (solvent) 
have  been  discussed  as  models  for  condensed  phase  systems.  ^Tnltial 
studies  of  these  models  have  been  focused  on  relatively  small  clusters 
with  few  solvent  molecules:  geometries,  energetics,  and  energy 

f  dynamics  have  all  been  explored  in  these  vdW  systems.. 

. . ^ 

In  a  companion  paper  (here  after  referred  to  as  l),  information 
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on  benzene  solvated  with  methane,  ethane,  and  propane  is  presented. 
Based  on  this  and  previous  work  from  our  laboratory,  the  nature  of 
solvation  and  the  solvation  process  has  been  explored.  The  solvent 
shift  in  the  benzene  clusters  is  understood  to  arise  primarily  from 
the  interaction  of  the  solvent  with  the  tt -cloud  of  the  aromatic  ring: 
hence,  the  strong  dependence  of  spectral  shift  on  geometry. 
Identification  of  cluster  geometry  and  energy  levels  has  also  led  to 
understanding  of  the  nucleation  processes  important  in  the  formation 


of  van  der  Waals  clusters.  Homogeneous  nucleation  can  occur  through 
the  vibrational  predissociation  (VP)  of  a  solvent  dimer  as  it  forms  a 
complex  with  the  aromatic  solute.  Inhomogeneous  nucleation  (the 
addition  of  two  or  more  solvent  molecules  at  a  time  to  a  solute)  can 
be  the  dominant  process  in  the  beam:  the  relative  amount  of 

•■d 


homogeneous  and  inhomogeneous  nucleation  is  proportional  to  the  ratio 
of  solvent-solvent  binding  energy  to  solvent-solute  binding  energy. 


□  □ 


The  present  paper  reports  the  results  of  a  study  of  the  vd¥ 
clusters  of  toluene  with  methane,  ethane,  and  propane.  Due  to  the 
reduction  of  the  solute  symmetry  compared  to  benzene,  the  geometry  and 
energetics  of  toluene  clusters  are  different  from  those  of  benzene 
clusters.  Comparison  of  the  benzene  and  toluene  results  for  similar 
beam  conditions  and  the  same  solvents  strengthens  and  generalizes  the 
conclusions  reached  previously. 

Both  large  (tol(s)^,  x>3)  and  small  (x<3)  clusters  have  been 
studied  for  the  toluene  solvent  systems.  The  changes  in  the  spectra 
of  these  clusters  as  larger  and  larger  clusters  are  probed  indicate 
possible  relationships  between  the  clusters  and  condensed  phase 
systems.  Comparisons  of  the  large  cluster  spectra  with  the  spectra  of 
cryogenic  solutions^  suggest  how  well  and  to  what  extent  the  jet 
generated  clusters  model  condensed  phase  systems. 

Computer  modeling  of  toluene-hydrocarbon  clusters  is  extremely 
helpful  in  supplementing  the  small  cluster  spectroscopic  data  in  order 
to  determine  cluster  geometry  and  binding  energy.  A  number  of  local 
minima  are  found  in  the  potential  surface  for  toluene  interacting  with 
these  hydrocarbon  solvents:  only  the  lowest  energy  computer  generated 
geometries  which  are  consistent  with  the  observed  spectra  are  taken  as 
acceptable  solute-solvent  configurations.  If  the  computer  generated 
potential  energy  surface  has  a  number  of  minima  with  relatively  low 
energy,  that  geometry  which  best  agrees  with  the  toluene  and 
benzene  spectra  is  accepted  as  the  appropriate  one. 
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The  toluene  cluster  data  will  be  discussed  in  two  parts. The  data 
for  small  clusters  will  be  analyzed  first  and  will  be  compared  to  the 
benzene  small  cluster  data.  The  modeling  of  condensed  phase  systems 
by  large  clusters  will  be  dealt  with  at  the  end  of  the  Discussion 
Section. 

II  EXPERIMENTAL  PROCEDURES 

The  apparatus  and  techniques  employed  in  this  work  are  similar  to 
those  described  in  I.  Mass  selective  absorption  (2-color  MS)  is  the 
only  spectroscopic  technique  used  in  this  study.  The  optimal  backing 
conditions  for  each  nozzle  are  found  experimentally  and  are  reported 
in  the  figure  captions. 

Three  different  nozzles  are  used  in  these  experiments.  The 
pulsed  nozzle  with  a  500  ym  orifice  produced  a  factor  of  '^15  larger 
intensity  than  the  50  urn  continuous  (CW)  nozzle  which  in  tuni 
produces  signals  3  to  4  times  more  intense  than  the  25Um  CW  nozzle. 
Greater  enhancements  in  the  signals  were  expected  for  the  larger 
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diameter  nozzles;  several  possible  loss  mechanism  may  contribute 
to  this  reduced  signal  improvement.  In  general,  however,  direct 
comparison  between  signals  from  the  various  nozzles  is  difficult 
because  the  optimum  backing  conditions  for  clustering  are  different 
for  each  nozzle. 

Lowering  the  backing  pressure  and  raising  the  pulsed  nozzle 
temperature  reduces  the  cooling  in  the  expansion  and  produces  hot 
clusters  which  evidence  spectral  hot  bands.  Although  the  (calculated) 
translational  temperature  of  the  expansion  at  the  most  extreme  warming 
conditions  (P^  =  10  psi,  T^  =  70°C)  is  still  quite  low  (  '''IK),® 
the  vibrational  temperature  of  the  clusters  is  ''<50  to  100K  based  on 
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the  observed  spectral  intensities.  Calculations  of  the  minimum  energy 
configurations  of  toluene  clusters  are  performed  as  described  in  I. 

The  coe.fficients  for  the  interactions  between  the  methyl  group  of 
toluene  and  the  solvents  are  taken  to  be  the  same  as  those  between 
solvents.  The  coefficients  are  reproduced  in  Table  I  for  convenience. 

Minimization  of  the  energy  of  a  given  cluster  is  accomplished  by 
moving  the  ligands  in  a  stepwise  fashion,  calculating  the  energy  at 
each  step,  and  comparing  the  two  most  recent  values  of  the  energy.  The 
procedure  will  terminate  each  time  a  local  minimum  in  the  energy  has 
been  found.  To  ensure  that  all  minima  are  found  and,  in  particular, 
that  the  lowest  minimum  is  located,  all  reasonable  starting  geometries 
must  be  employed.  Low  energy  and  agreement  with  the  spectroscopic 
data  are  the  two  criteria  used  to  settle  on  "correct"  configurations 
for  the  clusters. 

Ill  RESULTS 

Toluene-Methane 

The  Tol(CH^)^  2  absorption  spectra  have  been  published 

previously  and  are  reproduced  here  in  figure  1  and  Table  II  for 

convenience.  These  spectra  exhibit  sharp  origins  and  features  due  to  van 

der  Waals  (vdW)  stretching  and  bending  modes.  The  spectral  shift  for 

the  To1(CH^)2  Oq  is  twice  that  of  the  Tol(CH^)^  0^.  The 

additive  shift  indicates  addition  of  the  two  methane  molecules  to 

equivalent  sites  on  the  toluene.  The  stretching  mode  for  Tol(CH^)^ 
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has  a  considerable  anharmonicity  as  pointed  out  previously. 

To1(CH^)2  evidences  a  similar  progression  in  the  stretching  mode. 

The  To1(CH^)2  peak  at  -20  cm"^  is  probably  associated  with  a 
second  configuration  of  Tol(CH. )„. 


Larger  clusters  of  toluene  with  methane  produce  broad  spectra 

which  become  more  diffuse  and  structureless  as  the  clusters  increase 

in  size  (see  figure  2).  The  red  shift  of  the  most  prominent  feature 

in  the  Tol(CH^)^  spectrum  is  non-additive;  this  feature  is  actually 

less  red  shifted  than  the  TolCCH^)^  Oq.  Nonetheless,  intensity  is 

still  evident  in  the  Tol(CH^)^  spectrum  far  to  the  red  of  the 
0  5 

To1(CH^)2  Oq*  The  absorption  spectra  of  clusters  as  large  as 

To1(CH^)^q  have  been  observed.  The  spectra  become  very  broad  and 

unstructured  with  a  maximum  about  -50  cm~^  from  the  Tol  0^. 

The  binding  energy  of  these  clusters  can  be  estimated  by 

observing  the  presence  or  absence  of  higher  vibronic  bands  in  the 

2-color  MS  spectra.  If  the  first  photon  induces  a  transition  to  a 

level  which  results  in  vibrational  predissociation  of  the  cluster,  the 
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transition  will  not  be  observed  in  the  2-color  MS  spectrum.  For 
Tol(CH^)^  2  ■the  6bQ  transition  is  observed  but  the  1q 
transition  is  not.  This  sets  the  limits  on  the  binding  energy  at  554 £ 
Do1740  cm‘^  . 

The  binding  energy  is  the  same  within  these  limits  for 
Tol(CH^)^  2*  The  binding  energy  for  larger  clusters  is  difficult  to 
estimate  since  the  spectra  of  large  clusters  are  less  well  defined  and 
weak. 

The  computer  modeling  of  Tol(CH^)^  shows  a  minimum  energy 
configuration  in  which  the  methane  is  shifted  away  from  the  center  of 
the  aromatic  ring  toward  the  methyl  group  (see  figure  5,  Table  III). 

The  methane  should  interact  strongly  with  the  aromatic  it -cloud  of 
toluene  in  this  predicted  geometry.  Only  one  stable  configuration  is 
calculated  for  Tol(CH^)^:  this  agrees  with  the  observation  of  one 
origin  in  the  spectrum.  The  calculated  binding  energy  of  667  cm  ^  is 


consistent  with  the  limits  obtained  spectroscopically. 


Two  acceptable  geometries  are  found  in  the  computer  modeling  of 
To1(CH^)2.  The  configuration  labeled  I  in  figure  5  features  a 
methane  on  either  side  of  the  aromatic  ring,  forming  an  isotropic 
cluster.  The  second  configuration  predicted  by  the  computer  modeling 
(labeled  II  in  figure  5)  is  anisotropic  with  both  methane  molecules  on 
the  same  side  of  the  ring.  In  the  isotropic  geometry  both  of  the 
solvent  molecules  occupy  equivalent  sites  similar  to  the  site  occupied 
by  methane  in  Tol(CH^)^.  Consequently,  this  geometry  should  have 
the  larger  (additive)  spectral  shift  in  the  Tol(CH^)2  spectrum 
( figure  1 ) . 

Toluene-Ethane 

The  small  clusters  of  toluene  and  ethane  (Tol(C2Hg)^  2^ 
exhibit  sharp,  assignable  structure  (see  figure  4  and  Table  IV).  Tol 
^  long  progression  in  the  vdW  stretch,  as  well  as  two 
intense  bends.  The  Tol(C2Hg)2  spectrum  clearly  suggests  that  two 
geometries  exist  for  this  cluster.  One  geometry  generates  a  spectral 
shift  which  is  twice  that  found  for  Tol(C2Hg)^.  The  vdW 
vibrational  features  of  this  additive  shift  configuration  consist  of 
long  progressions  in  the  stretching  mode.  The  second  geometry  is 
characterized  by  a  small  (-10  cra”^ )  solvent  shift  and  a  short 
progression  of  doublets  after  the  origin  peak.  The  absorption  spectra 
of  larger  Tol(C2H^)^  clusters  become  broad  with  very  little 
well-defined  structure.  The  maximum  in  a  typical  absorption  spectrum 
occurs  at  about  50  cm”^  to  the  low  energy  side  of  the  Tol  Oq, 

although  significant  measurable  intensity  extends  to  roughly  -200 

-1 

cm 

The  effect  of  warming  the  Tol(C2Hg)^  clusters  is  similar  to 
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that  of  warming  the  Tol(CH^)^  nluaters.  The  hot  hands  produce  a 
broad  background  somewhat  to  the  low  energy  side  of  the  cold  bands.  As 
judged  by  the  decrease  in  the  signal  intensities,  continued  warming  (P 
<  20  psi,  T^  '^70°C)  begins  to  destroy  the  vdW  complexes. 

Vibrational  predissociation  of  TolCC^Hg)^  2  dusters  occurs 
at  the  12^  toluene  excitation  but  not  at  the  binding  energy  of 

these  clusters  therefore  lies  between  740  and  933  cm  ^ .  Both 
geometries  of  Tol(C2Hg)2  seem  to  have  the  same  binding  energy 
within  this  experimentally  defined  range. 

The  computer  modeling  of  Tol(C2Hg)^  shows  two  minima  in  the 
potential  surface  above  toluene.  Since  only  one  geometry  is  obvious 
experimentally  only  the  lower  energy  computer  generated  configuration 
is  presented  in  figure  5  and  Table  V.  The  less  stable  geometry  (  '\;180 
cra~^  higher  in  energy)  is  probably  not  observed  spectroscopically, 
although  one  of  the  peaks  assigned  to  a  vdW  bending  mode  is  possibly 
the  origin  of  a  second  configuration.  The  binding  energy  predicted  by 
the  computer  modeling  is  within  the  spectroscopically  determined 
limits. 

For  Tol(C2Hg)2  several  possible  geometries  are  calculated. 

Si:ice  the  spectrum  clearly  shows  only  two  geometries,  the  most 
probable  configurations  must  be  chosen  from  this  group.  One  of  the 
geometries  must  be  an  isotropic  configuration  with  an  ethane  on  either 
side  of  the  ring  and  the  other  configuration  must  involve  both  ethane 
molecules  on  the  same  side  of  the  ring.  Configurations  consistent 
with  the  spectrum  are  shown  in  figure  5. 


Toluene-Propane 

The  spectra  of  the  small  clusters  of  toluene  and  propane 
(To1(C_Ho)<  o  i)  also  consist  of  well  resolved  features  (figure  6 

y  O  1  ^ 

and  Table  Vl).  The  To1(C^Hq);|  cluster  spectrum  strongly  suggests 
two  configurations  exist  for  this  cluster.  The  geometry  with  a 
relatively  large  spectral  shift  shows  a  progression  in  the  vdW  stretch 
as  well  as  several  vdW  bends  to  the  high  energy  side  of  its  0^ 
transition.  The  second  configuration  gives  rise  to  a  small  spectral 
shift  of  the  origin:  the  origin  is  followed  by  a  progression  of 
doublets.  This  pattern  is  quite  reminiscent  of  that  found  for  the 
second  (small  shift)  configuration  of  Tol(C2Hg)2* 

The  vdw  cluster  spectrum  of  To1(C^Hq)2  has  a  complicated 
structure  indicating  three  or  four  possible  different  configurations 
for  this  mass  species.  These  geometries  are  generally  distinguishable 
baaed  on  relative  intensities,  spectral  shifts,  qualitative 
vibrational  analysis,  and  comparison  with  the  Tol(CH^)^  2> 

Tol(C2Hg)^  2*  Tol(C^Hg)^  data. 

The  Tol(C^Hg)^  spectrum  is  actually  less  congested  than  the 
Tol(C^Hg)2  spectrum,  presumably  because  some  of  the  possible 
configurations  are  not  observed  in  the  beam.  Two  sets  of  peaks, 
corresponding  to  two  different  geometries  of  Tol(C^Hg)^,  are 
present  in  the  spectrum  as  shown  in  figure  6.  One  set  of  peaks 
features  a  regular  progression  in  an  apparently  single  vdW  mode  and 
the  other  features  a  short  progression  of  doublets. 

Warming  the  toluene-propane  clusters  destroys  the  well-resolved 
structure  in  the  Tol(C^Hg)2  ^  spectra  (figure  7);  as  in  the  other 
systems,  hot  bands  of  the  Tol(C^Hg)^  clusters  form  an  unstructured, 
broad  background  to  the  red  of  the  cold  spectrum. 

"’he  large,  clusters  of  toluene  and  propane  produce  spectra 
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qualitatively  similar  to  those  of  the  other  systems  (figure  8). 
Although  some  structure  exists  in  the  ToI(C^Hq)^  spectrum,  the 
Tol(C^Hg)g  spectrum  is  largely  unstructured  with  an  intensity 
maximum  around  -60  cm”^  relative  to  the  Tol  Og. 

The  binding  energies  of  the  Tol(C^Hg)^  2  3  clusters  are  all 
greater  than  933  cm~^  as  all  three  clusters  show  clear  spectra  in 
the  1 2g  region.  No  toluene  features  above  12q  are  either  intense 
enough  or  free  enough  from  interferences  to  allow  further  studies; 
therefore,  no  upper  limit  to  the  binding  energy  can  be  quoted  at  this 
time . 

Computer  modeling  predicts  several  stable  geometries  for 
Tol(C^Hg)^.  As  noted  earlier,  the  Tol(C^Hg)^  spectrum  indicates 
that  two  configurations  with  very  different  spectral  shifts  exist  (see 
figure  6).  Note  the  similarity  in  the  spectral  shifts  of 
configuration  II  of  Tol(C^Hg)^  and  configuration  II  of 
Tol(C2Hg)2  (figure  4).  This  similarity  in  spectral  shift  indicates 
a  similarity  in  the  overlap  of  the  solvent  with  the  aromatic  ring  of 
toluene.  Figure  9  and  Table  VTI  present  the  lowest  energy  geometries 
which  should  give  rise  to  spectral  shifts  consistent  with  the 
Tol(C^Hg)^  spectrum  base  on  the  correspondence  between 
solvent-aromatic  ring  overlap  and  spectral  shift  (l).  Other  less 
stable  geometries  may  be  producing  some  of  the  observed  peaks  in  the 
spectrum,  but  further  assignments  cannot  be  made  with  the  existing 
data. 

Similar  considerations  are  used  to  choose  those  geometries  of 
Tol(C^Hg)2  consistent  with  the  spectrum.  These  minimum 

energy  geometries  are  pictured  in  figure  10  and  tabulated  in  Table 
VIII.  Computer  modeling  of  To1(C,Hq)t.  was  not  done. 


IV  DISCUSSION 


The  spectral  shift  of  a  cluster  depends  primarily  on  the 
polarizability  of  the  solvent  along  the  intermolecular  bond  and  on  the 
change  in  polarizability  of  the  solute  upon  excitation.  One  would 
expect  that  the  change  in  polarizability  upon  excitation  is  larger  for 
toluene  than  benzene  and  thus,  that  the  cluster  shifts  are  larger  for 
toluene  than  benzene  clusters.  However,  most  of  the  change  in 
polarizability  upon  (  tr-Tr*  )  excitation  occurs  in  the  aromatic  ir-system; 
moving  the  solvent  away  from  the  ir -system  must  significantly  reduce 
the  cluster  shift.  Thus,  two  competing  effects,  increased 
polarizability  change  and  reduced  interaction  with  the  aromatic 
^-system,  influence  the  relative  spectral  shifts  of  benzene  and  toluene 
clusters  with  small  hydrocarbons.  In  the  ensuing  discussion  we  will 
treat  large  and  small  clusters  separately  and  discuss  shifts, 
geometries,  and  nucleation  processes. 

A.  Small  Clusters. 

Methane 

The  computer  modeling  of  the  Tol(CH^)^  cluster  shows  that  the 
methane  is  displaced  0.92lX  towards  the  methyl  group  of  toluene  (see 
figure  3,  Table  III).  This  calculated  geometry  presents  a  rational 
for  the  essential  (unexpected)  equivalence  of  the  cluster  spectral 
shifts  for  Tol(CH^)^  and  Ben(CH^)^,  even  though  toluene  has  a 
(presumed)  larger  change  in  polarizability  upon  excitation  than  does 
benzene.  Thus,  the  above  mentioned  two  competing  effects  tend  to 
cancel  one  another  for  Tol(CH^)^. 

These  basic  ideas  also  help  to  rationalize  the  similarities  and 


differences  in  the  TolCCH^)^  and  BenCCH^)^  spectra.  Spectra  of 
the  isotropic  clusters  (a  solvent  on  each  side  of  the  aromatic  ring) 
of  To1(CH^)2  and  BenCCH^)^  are  nearly  identical,  as  would  be 
surmised  based  on  the  above  arguments  for  the  Tol(CH^)^  clusters. 
However,  anisotropic  clusters  of  toluene  and  benzene  with  two  methanes 
exhibit  significantly  different  shifts.  Computer  modeling  of  the 
anisotropic  cluster  of  Ben(CH^)2  suggests  that  only  one  methane 
moiety  interacts  strongly  with  the  tt -system  (l),  thus  yielding  a 
spectral  shift  similar  to  that  observed  for  Ben(CH^)^.  On  the  other 
hand,  computer  modeling  of  the  anisotropic  cluster  of  Tol(CH^)2 
shows  that  both  methane  molecules  are  considerably  displaced  from  the 
center  of  the  aromatic  ring  (1.09^  and  2. ’37%)-  The  observed  small 
spectral  shift  for  TolCCH^)^  (-21cm”^)  is  thereby  consistent  with 
the  previously  determined  (I)  relationship  between  spectral  shift  and 
geometry. 

The  ratios  of  intensities  of  the  isotropic  to  anisotropic  cluster 
peaks  are  different  for  TolCCH^)^  and  BenCCH^)^:  ^iso^^aniso 

[To1(CH.  )_,0?]  '^5  (figure  1),  whereas  I  .  /I 

[Ben(CH^)2,  6q]  M .  As  pointed  out  in  I,  isotropic  clusters  can 
only  be  produced  by  homogeneous  nucleation.  Most  homogeneously 
nucleated  clusters  are  argued  to  form  through  the  collision  of  a 
solvent  dimer  with  the  solute;  the  binding  energy  of  the  solute  with 
one  of  the  solvent  molecules  is  dissipated,  in  part,  by  VP  of  the 
solvent  dimer  bond.  Apparently,  the  larger  solute-solvent  binding 
energy  of  the  toluene  clusters  (as  predicted  by  the  computer  modeling) 
facilitates  breaking  of  the  solvent  dimer  bond  and  allows  homogeneous 


nucleation  to  dominate  the  formation  of  Tol(CH^)2  clusters. 

Nucleation  processes  will  be  treated  in  more  detail  at  the  end  of  this 
section. 

Ethane 

Figure  5  depicts  the  computer  calculated  geometry  for 
Tol(C2Hg)^:  the  ethane  molecule  is  displaced  toward  the  methyl 
group  of  toluene  and  inclined  with  respect  to  the  plane  of  the  ring. 
Assuming  this  geometry  to  be  correct,  it  results  in  a  spectral  shift 
smaller  by  roughly  5  cm  ^  than  that  found  for  the  comparable 
Ben(C2Hg)^  configuration.  The  methyl  group  of  toluene,  according  to 
the  potential  calculations,  draws  the  ethane  solvent  molecule  away 
from  the  ring  center,  thereby  reducing  the  toluene  spectral  shift 
relative  to  the  comparable  benzene  cluster. 

An  isotropic  Tol(C2Hg)2  cluster  is  clearly  indicated  in  the 
observed  spectrum  based  on  an  origin  shift  twice  that  of  the 
Tol(C2Hg)^  origin  (figure  4).  The  intensity  of  features  due  to 
isotropic  clusters  relative  to  those  due  to  anisotropic  clusters  is 
much  less  than  seen  in  the  Tol(CH^)^  spectrum.  The  ratio  is 
consistent  with  the  trend  report  in  I;  homogeneous  nucleation  can  be 
less  probable  with  solvents  having  larger  binding  energies  than  with 
solvents  having  smaller  binding  energies.  However,  more  homogeneous 
nucleation  is  evident  for  the  toluene-ethane  system  than  for  the 
benzene-ethane  system;  this  same  trend  was  found  for  the 
solute-methane  clusters  and  is  related  to  the  strength  of  the 
solvent-solute  bond. 
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•  ^ 

^  •  I 


/-  -V 


Propane 


Arguments  put  forth  above  concerning  spectral  shift, 
solvent-solute  geometry,  binding  energies,  etc.,  can  also  be  applied 
for  the  somewhat  more  complicated  toluene-propane  clusters.  The 
spectral  shifts  of  the  two  configurations  determined  by  computer 
modeling  of  Tol(C^Hg)^  (figure  9)  can  thereby  be  qualitatively 
predicted.  Configuration  I,  having  a  significant  interaction  of  a 
terminal  methyl  group  of  propane  with  the  aromatic  it-system  of 
toluene,  should  have  a  spectral  shift  greater  than  that  of 
Tol(CH^)^ .  On  the  other  hand,  configuration  II  (figure  9)  would  be 
expected  to  have  a  small  spectral  shift  since  the  propane  molecule  is 
close  to  the  methyl  group  of  toluene.  Origins  with  the  expected 
shifts  are  indeed  observed  in  the  spectrum  of  Tol(C^Hg)^  as 
depicted  in  figure  6. 

The  Tol(C^Hg)2  spectrum  (figure  6)  has  substantial  intensity 
tentively  assigned  to  isotropic  clusters  (see  origin  peaks  labled  I, 
II,  III  in  Table  VI  and  figure  10:  although  it  is  also  possible  that 
configurations  II  and  III  may  be  anisotropic.  This  observation 
apparently  does  not  fit  the  trend  established  above  if  the  computer 
predicted  values  of  the  solvent  dimer  and  solute-solvent  binding 
energies  are  used.  However,  computer  modeling  of  the  propane  dimer 
suggests  that  several  geometries  are  stable,  some  of  which  have 
relatively  small  binding  energies.  If  a  reasonable  fraction  of  the 
propane  dimers  in  the  beam  are  hot  or  are  in  geometries  with  small 
binding  energies  (i.e.  ca.  500  cm~^),  the  relative  number  of 
isotropic  clusters  in  the  beam  would  be  larger  than  otherwise 
expected. 
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Nucleation 


The  nucleation  process  in  the  jet  expansion  can  be  better 
understood  by  employing  a  simple  kinetic  scheme.  We  assume,  as  has 
been  previously  suggested  (l),  that  small  vdW  clusters  are  formed  by 
solute  collisions  with  solvent  dimers  (or  possibly  larger  species). 
While  this  may  not  be  the  only  cluster  formation  mechanism,  it  is 
consistent  with  all  benzene  and  toluene  cluster  data  and  can  be 
justified  by  reasonability  arguments.  A  solvent  dimer  colliding  with 
the  solute  produces  an  "activated  complex"  in  much  the  same  fashion  as 
is  typically  considered  in  a  bimolecular  reaction.  This  complex  can 
then  decay  through  three  different  pathways:  I )  vibrational 
predissociation  of  the  solvent  dimer  to  form  a  homogeneously  nucleated 
solute  (solvent)^  vdW  cluster;  2)  collision  with  a  third  body  (He, 
typically)  to  generate  a  inhomogeneous ly  nucleated  solute  (solvent)^ 
vdW  cluster;  or  3)  decay  back  to  the  reactants  (solvent)^  and  solute. 

Data  presented  in  the  Results  Section  indicate  that  details  of 
the  solvent  dimer  binding  energy  and  solute-solvent  binding  energy  are 
important  for  the  relative  concentrations  of  homogeneously  and 

inhomogeneous ly  nucleated  products.  Thus,  finer  detail  in  the  overall 
kinetic  scheme  should  be  considered.  The  mechanism  suggested  for 
the  reaction  is  represented  as  follows: 


I  II 


TT 


0—0  + 


in  which indicates  a  vdW  bond  with  significant  vibrational 
excitation,  0  indicates  a  solvent  molecule,  M  stands  for  an  inert 
third-body  collision  partner  and  - —  represents  the  solute.  Assuming 
d[ll]/dt=0,  and  ignoring  the  first  step  in  the  reaction  chain,  the 
ratio  of  concentrations  [iv]  to  [ill]  is  given  by 


CIV3  (k^4k2-+k3-Hc^) 

—  ■  ig  q 


(1) 


The  amount  of  inhomogeneous  nucleation  relative  to  homogeneous 
nucleation  is  proportional  to  this  ratio.  Two  reasonable 
approximations  can  be  made  for  this  expression.  First,  <k^ ,k2 
as  k^  is  at  most  sec”^  in  the  backing  region  (and 

certainly  much  leas  in  the  beam)  while  k^  ,k2  ''^10  sec  .  If 
this  were  not  the  case,  inhomogeneously  nucleated  clusters  would 
always  dominate  the  intensity.  Second,  k^s  k2  because  the  density 
of  states  is  high  for  both  species  at  the  given  energies  (500-900 
cm~^ )  and  energy  is  not  expected  to  concentrate  in  one  bond  or  the 
other  over  many  vibrational  periods.  With  these  approximations 
equation  (1)  becomes 

CIV3  k^ 

[III]  ^ 


k^  depends  only  on  the  collision  rate  of  third  bodies  with  the 
species  I  and  II  and  should  be  roughly  the  same  for  all  the  systems 
studied  here.  The  difference  in  the  relative  numbers  of  isotropic 
species  for  the  various  systems  must  then  be  associated  with  the 
changes  in  k^  from  cluster  to  cluster. 
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The  transformation  of  species  II  to  species  III  can  be  viewed 
as  a  reaction  with  a  potential  barrier  equal  to  the  binding  energy  of 
the  solvent  dimer.  The  rate  constant  can  be  written  as 


kj  •V  exp  [ 


in  which  is  the  solvent  dimer  binding  energy  and  is  the 

energy  in  the  solvent  dimer  bond  for  species  II.  If  equilibrium  is 
established  (k^,k2>>k^)  before  VP  occurs  and  the  excess  energy  in 
the  intermediate  is  statistically  distributed  between  the 
solvent-solute  and  solvent-solvent  bonds,  then  Ej^j^ 
which  BEgj^  is  the  solute-solvent  binding  energy.  This  latter 
quanity  takes  on  the  role  of  kT  for  the  reaction  II'*'  III.  Thus  k^ 
becomes  a  function  of  both  the  solvent-solvent  and  solute-solvent 
binding  energies.  Table  IX  shows  that  the  relative  number  of 
isotropic  clusters  in  the  system  scales  approximately  with  this  form 
for  k^.  Of  course,  the  detailed  behavior  of  the  ratio  [lll]/[ll] 
depends  on  the  appropriateness  of  the  various  approximations  for  each 
system. 

B.  Large  Clusters 

The  intensity  maximum  in  the  spectrum  of  clusters 

with  more  than  three  solvent  molecules  is  much  less  red  shifted  with 

respect  to  the  free  molecule  than  the  intensity  maximum  in 

5 

the  comparable  cryogenic  solution.  This  observation  is  directly 


attributable  to  the  conditions  of  the  experiment  and  the  nature  of  the 
nucleation  processes  in  the  expansion.  Inhomogeneous  nucleation 
dominates  the  formation  of  large  clusters,  especially  for  clusters 


involving  better  solvents  such  as  propane.  Since  inhomogeneous 
nucleation  produces  more  anisotropic  clusters,  most  of  the  large 
clusters  formed  in  the  jet  will  be  anisotropic  and  will  exhibit  small 
spectral  shifts;  larger  spectral  shifts  are  found  for  isotropic 
(solution-like)  clusters. 

The  spectral  shifts  in  the  jet  generated  clusters  have  been 
determined  to  be  largely  due  to  interaction  between  the  solvent  and  the 
aromatic  ir-system  of  the  solute.  If  the  clusters  are  to  serve  as 
potential  models  of  condensed  phase  systems,  the  difference  between 
the  principal  gas  to  cluster  and  gas  to  liquid  shifts  must  be 
explored.  [Note  nonetheless,  that  some  cluster  intensity  is  found  at 
the  same  spectral  shift  as  that  observed  for  the  analogous 
solution.  ’  ’  J  While  in  general  the  shifts  for  clusters  are 
saturated  with  a  solvent  molecule  on  each  side  of  the  aromatic  ring, 
larger  shifts  might  arise  through  two  mechanisms:  forced  crowding  of 
the  solvent  around  the  solute  so  that  more  molecules  effectively 
interact  with  the  tr -system  of  the  solute,  and  interaction  of  the 
solvent  with  the  solute  in  the  plane  of  the  ring.  These  latter 
positions  are  not  local  minima  in  the  potential  surface  of  small 
clusters.  However,  the  cage  structure  of  the  liquid  would  force 
solvent  molecules  to  interact  with  the  solute  in  the  ring  plane.  This 
in-plane  interaction  may  also  contribute  to  the  spectral  shift  of 
larger  vdW  clusters  and,  in  particular,  could  be  responsible  for  the 
weak  intensity  at  approximately  the  value  of  the  solution  spectral 
shift  reported  for  these  systems  previously.^  Since  we  have  no  direct 
information  on  the  in-plane  interactions,  the  relative  contribution  to 
the  shift  from  these  two  possible  sources  is  difficult  to  assess. 


V. 


CONCLUSIONS 


Two-color  time  of  flight  mass  spectroscopy  and  simple  potential 
calculations  have  been  employed  to  determine  cluster  binding  energy, 
spectral  shift,  geometry  and  nucleation  processes  for  toluene  solvated 
by  methane,  ethane  and  propane.  The  major  findings  of  this  effort  can 
be  summarized  as  follows: 

1 .  Cluster  binding  energy  is  related  to  the  polarizability  of 
the  two  species  in  the  cluster  and  is  apparently  not  a  sensitive 
function  of  solute-solvent  geometry; 

2.  The  observed  spectral  shift  for  a  cluster  is  due  mostly  to 
the  solvent  interaction  with  the  solute  aromatic  -system  and 
therefore  is  a  sensitive  function  of  cluster  geometry; 

3.  The  nucleation  of  clusters  under  the  experimental  beam 
conditions  is  suggested  to  occur  through  the  interaction  of  solvent 
dimers  or  larger  aggregates  and  the  solute  monomer  -  whether  isotropic 
clusters  (homogeneous  nucleation)  or  anisotropic  clusters 
(inhomogeneous  nucleation)  arise  depends  on  the  relative  size  of  the 
solvent  dimer  binding  energy  and  the  solute-solvent  binding  energy; 

4*  The  benzene  and  toluene  cluster  systems  appear  to  behave  in 
a  consistent  and  predictable  fashion  -  the  data  analysis  for  both 
systems  has  been  carried  out  such  that  trends  are  transferable  from 
one  system  of  clusters  to  another;  and 

5.  The  exponential-6,  atom-atom  potentials  employed  in  the 
computer  modeling  of  cluster  geometry  and  binding  energy,  which  have 
been  quite  successfully  used  over  the  years  in  treating  condensed 
phases,  are  apparently  reasonable  potentials  for  the  modeling  of 
cluster. 
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TABLE  I 


Parameters  for  the  energy  esqnression  in  the  conpater  oodeling  of  the  vdW  clusters 
of  toluene. 


^ij  *  -  /^ij  +  Bij  exp  [  -  Cy 


-1 

p/  cm  V 

c  (8  ‘h 

role  ' 

'  role  ^ 

Aliphatic-Aliphatic® 

C-C 

1310% 

8890353 

3.A21 

C-H 

A7830 

7562708 

3.9A 

H-H 

15028 

6390528 

A.6A3 

A1  ipha  t  ic  -  ArooH  tic^ 

C-C 

1626A5 

16520286 

3.5105 

C-H 

4786A 

6A3A085 

3.805 

H-H 

13071 

2A97533 

A.  1915 

Aroma  tic -Aromatic^ 

C-C 

201786 

30698A27 

3.6 

C-H 

A7897 

5A73891 

3.67 

H-H 

11368 

976084 

3.7A 

From  the  metlane-inethane  parameters  of  ref.  11  in  paper  I. 

The  parameters  are  fovnd  by  combining  the  benzene-benzene  and  methane-methane 
parameters  as  follows: 


Sm  W 

From  the  benzene-benzene  parameters  of  ref.  10  in  paper  I. 


TABLE  II 


I 

ir 


I? 

kil 


Tw-color  TOFMS  of  Tol(CH^)  and  Tol(CH^)2  In  the  o”  region.  Only  the  prominent  peaks 
are  tabulated  and  tentative  assignments  are  given  (see  figure  1). 


Tol(CH^) 


37434.4 
37451.8 
37458.7 
37466.2 

37488.5 


To1(CH^)2 


Tentative  assignments  are  as  follows:  V  stands  for  the  Tol-CH,  stretch 
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and  A  and  B  are  different  vdW  bends.  I  and  II  represent  different 
configurations. 


TABLE  111 


The  atomic  positions  for  the  minimLin  energy  configurations  of  Tol  (CH^)  j^and  Tol(CH^)2 

(see  figure  3).  The  origin  is  at  the  center  of  the  aromatic  ring  of  toluene.  The 
+  X  axis  passes  through  the  methyl  group,  the  Y  axis  is  in  the  plane  of  the  ring 
and  the  Z  axis  is  out  of  the  plane  of  the  ring- 


IP  - 

*  -  • 

•  -*• 

Atom 

x(X) 

Y(X) 

z(X) 

y/- 

1 

Toluene 

Cl 

2.925 

0 

0 

C2 

-1.395 

0 

0 

'  k  *• 

C3 

-  .6975 

1.208 

0 

l-kk. 

• 

C4 

.6975 

1.208 

0 

C5 

1.395 

0 

0 

C6 

.6975 

-1.208 

0 

C7 

-  .6975 

-1.208 

0 

HI 

3.29 

0 

-1.038 

H2 

3.29 

.8985 

.518 

H3 

3,29 

-  ,8985 

.518 

1 

HA 

-2.A79 

0 

0 

t: 

H5 

-1.2A 

2.147 

0 

m 

H6 

1.2A 

2.147 

0 

H7 

1.2A 

-2.147 

0 

a 

H8 

-1.2A 

-2.147 

0 

Tol(CH^)^ 

Cl 

.921 

0 

3.384 

HI 

-  .104 

0 

2.985 

H2 

1,451 

-  .8985 

3.036 

•: 

H3 

1.A51 

-  .8985 

3.036 

HA 

.884 

0 

4.483 

»  •'  , 

To1(CH^)2 

Configuration 

I 

Ligand  1 

Cl 

.922 

0 

3.381 

HI 

-  .103 

0 

2.982 

H2 

1.452 

-  .8985 

3.033 

H3 

1.452 

.8985 

3.033 

HA 

.886 

0 

4.48 

I---. 

V 

F*  ^  * 

Ligand  2 

Cl 

.767 

0 

-3.391 

HI 

-  .255 

0 

-2.983 

»'  ^  -* 

H2 

1,300 

-  .8985 

-3.048 

• .  • 

H3 

1,300 

.8985 

-3.048 

* « • 

HA 

.721 

0 

-4.49 

Energy 
(cm  ) 


-  667 


-1348 


‘•VAv'.'A-/ 


TABLE  III 
(continued) 


Atom 

X(X) 

Y(X) 

Z(X) 

Energy 

(cm"^) 

Configuration  11 

-1360 

Ligand  1 

Cl 

2.372 

.003 

3.433 

HI 

1.977 

-  .8968 

2.938 

H2 

3.470 

.005 

3.371 

H3 

1.973 

.900 

2.938 

H4 

2.066 

.003 

4.490 

Ligand  2 

Cl 

-1.089 

-  .003 

3.454 

HI 

-2.172 

-  .005 

3.260 

H2 

-  .634 

-  .901 

3.011 

H3 

-  .638 

.896 

3.011 

HA 

-  .913 

-  .003 

4.540 

s  > 


TABLE  IV 


Two-color  TOFMS  of  Tol(C2H^)j  and  Tol(C2H^)2*  Only  the  prominent  peaks  are  tabulated 
and  tentative  assignments  are  given  <see  figure  4). 


Species 


Energy 


(vac  cm 


Energy  origin  relative  Energy  relative 


Assignment 


to  Tol  Oq  (cm 


to  cluster  0^  (cm  S 


Tol(C  H  ) 

37435  -43 

.0 

37450 

15 

37454 

19 

37458 

23 

< 

37460 

25 

*0 

37475 

40 

37485 

50 

37486 

51 

'>1 

Tol(C-H,).  I 

*^Tentatlve  assignments  are  as  follows:  V  stands  for  the  fundamental  vdW  stretch 
A  and  B  are  different  vdW  bends.  I  and  II  represent  different  geometries.  In 


Tol(C,H,)-  the  vdW  modes  may  involve  solvent-solvent  motions  as  well  as 
2  o  2 


solvent-solute  motions. 


TABLE  V 


The  atomic  positions  for  the  minimum  energy  configurations  of  Tol(C_H,),  and 

Z  O  I 

Tol(C«H,)_  (see  figure  5).  Axes  and  solute  coordinates  are  the  same  as  in  Table  111 


Tol<CjHj)j 

Ligand 


Ligand  1 


Ligand  2 


Configuration  II 
Ligand  1  C 


Cl 

.166 

-  .001 

3.742 

C2 

1.672 

-  .001 

3.478 

HI 

-  .374 

0 

2.783 

H2 

-  .104 

.897 

4.317 

H3 

-  .105 

-  .900 

4.315 

H4 

2.212 

-  .002 

4.437 

H5 

1.943 

.898 

2.905 

H6 

1.942 

-  .899 

2.903 

nfiguration  I 

Cl 

.017 

.041 

3.877 

C2 

1.515 

.185 

3.599 

HI 

-  .524 

-  .051 

2.924 

H2 

.  -  .340 

.929 

4.418 

H3 

-  .156 

-  .857 

4.487 

H4 

2.056 

.277 

4.553 

H5 

1.688 

1.083 

2.989 

H6 

1.574 

-  .734 

3.061 

Cl 

.302 

0 

-3.845 

C2 

1.802 

0 

-3.541 

HI 

-  .262 

.001 

-2.900 

H2 

.048 

.898 

-4.426 

H3 

.047 

-  .899 

-4.425 

H4 

2.366 

-  .001 

-4.486 

H5 

2.057 

.899 

-2.961 

H6 

2.056 

-  .898 

-2.960 

II 

Cl 

-  .628 

1.561 

3.933 

C2 

-1.263 

.311 

3.555 

HI 

-1.355 

2.254 

4.496 

H2 

-  .336 

2.187 

3.018 

H3 

.261 

1.471 

4.554 

H4 

-1.554 

-  .226 

4.469 

H5 

-2.151 

.490 

2.933 

H6 

-  .535 

-  .292 

2.991 

-  967 


-1663 


-1834 


TABLE  V 
(continued) 


X(X)  Y(X) 


Atom 


TABLE  VI 


Two-color  TOFMS  of  TolCC^Hg)^  ,  Tol(C2Hg)2  and  Tol(CgHg)g.  Only  the  prominent  peaks 
are  tabulated  and  tentative  assignments  are  given  (see  figure  6) . 


Energy 
(cm  S 


Energy  origin  relative  Energy  relative 
to  Tol  0^  (cm  to  cluster  0^  (cm 


Species 


Assignment 


TABLE  VI 
(continued) 


Species 

Energy 

(cm  ) 

Energy  origin  relative 
to  Tol  Oq  (cm 

Energy  relative 
to  cluster  0^  (cm 

Assignment^ 

11 

37A36 

-41 

0 

0? 

37A56 

14 

vj 

III 

37457 

-20 

0 

“S 

37479 

22 

vj 

IV 

37472 

-  5 

0 

37489 

16 

37492 

20 

TolCC^Hg)^  I 

37421 

-57 

0 

37434 

13 

< 

37446 

25 

II 

37454 

-24 

0 

< 

37469 

15 

4 

37477 

23 

4 

TABLE  VII 

The  atomic  positions  for  the  minimum  energy  configurations  of  ToI(C-Hp), 

3  o  1 

(see  figure  9).  The  origin  and  the  coordinates  of  the  solute  are  the  same  as  in 
Table  III. 


Atom  X(X)  Y(X)  Z(X)  Energy 

(cm”^) 


TABLE  VI 11 


The  atomic  positions  for  the  minimum  energy  geometry  of  configuration  IV  of 
ToKC^Hg)^  (see  figure  10).  The  origin  and  coordinates  of  the  solute  are  the 

same  as  in  Table  III. 


1 

1.9 

2 

.9 

3 

1.6 

1 

l.A 

2 

2.4 

3 

2.4 

o 


12 

1.278 

3.661 

4 

.155 

3.643 

i7 

-1.198 

3.633 

10 

2.251 

3.668 

10 

1.131 

4.435 

:9 

1.147 

2.888 

:0 

.230 

4.414 

:9 

.246 

2.866 

10 

-2.006 

3.620 

6 

-1.183 

4.410 

15 

-1.167 

2.862 

10 

.898 

3.577 

lO 

-  .490 

3.439 

i3 

-1.560 

3.585 

12 

1.665 

3.472 

.987 

4.566 

i3 

1.034 

2.792 

13 

-  .580 

2.450 

18 

-  .627 

4.224 

.6 

-2.558 

3.486 

15 

-1.423 

2.800 

-1.471 

4.574 

FIGURE  CAPTIONS 


FIGURE  1 

Two-color  TOPMS  of  ToUCH^)^  (upper  trace)  and  Tol(CH^)2  in  tne  Oq  region. 
The  energy  scale  is  relative  to  the  Tol  0^  (37477»5  cm”^).  Assignments  are 
given  in  Table  II.  This  figure  is  reproduced  from  reference  5* 

FIGURE  2 

Two-color  TOPMS  of  toluene-methane  large  clusters  in  the  0^  region.  The 

energy  scale  is  relative  to  the  Tol  0^  (37477.5  cm“^).  The  concentration  of 

methane  in  the  beam  is  ]%  for  Tol(CH.)_  and  2.5%  for  Tol(CH.),  and 

4  3  4  4 

Tol(CH^)^.  The  pulsed  nozzle  is  cooled  to  4°C  and  =  100  psia. 

FIGURE  3 

Minimiim  energy  configurations  of  Tol(CH^)^  and  Tol(CH^)2»  Atomic 
positions  are  given  in  Table  III.  Note  that  both  isotropic  (l)  and 
anisotropic  (II)  clusters  of  Tol(CH^)2  are  predicted  in  the  computer 
modeling. 

FIGURE  4 

Two-color  TOFMS  of  toluene-ethane  clusters  in  the  0^  region.  The  energy 
scale  IS  relative  to  the  Tol  0^  (37477-5  cm“^).  The  concentration  of  ethane 
in  the  beam  is  0.5%  for  Tol(C„H, ).  and  Tol(C_H, )„  and  2.0%  for  Tol(C„H, ),. 

4101  0  ^  d  O  J 

The  backing  conditions  for  the  pulsed  nozzle  are  P  =100  psia  and  T  =  4°C. 

o  o 

Data  for  Tol(C^H,).  and  Tol(C„H, )„  are  presented  in  Table  IV. 

^  O  1  ^  o  ^ 


FIGURE 


Minimum  energy  configurations  of  Toi(C„H^h  and  Tol(,C„H,  Atomic 

d  o  \  d  o  d 

positions  are  given  in  Table  V.  Tol(.C2H^)2  also  shows  isotropic  and 
anisotropic  clusters. 

FIGURE  b 

Two-color  TOFMS  of  toluene-propane  small  clusters  in  the  0^  region.  The 
energy  scale  is  relative  to  the  Tol  0^  (.37477. b  cm~^  j .  The  concentration  of 
propane  in  the  beam  is  \'jt>  for  all  tnree  spectra.  The  pulsed  nozzle  conditions 
are  -  100  psia  and  T^  =  4^C.  Data  for  these  spectra  are  presented  in 
Table  VI. 

FIGURE  7 

Two-color  TOFMS  of  Toi(,G,Huh  and  TolCC,HoJ„  under  hot  nozzle  conditions. 

3  o  1  3  o 

The  energy  scale  is  relative  to  the  Tol  0^  (37477.5  cm~^}.  The  concentration 
of  propane  is  1$,  =  40  psia  and  T^  70°C.  Note  the  dramatic  increase 

in  the  background  relative  to  the  spectra  in  Figure  6.  This  broad  background 
IS  due  to  hot  bands.  The  negative  peak  at  0  cm”^  in  the  Toi(C 

3  o  1 

spectrum  is  due  to  detector  overload  at  Tol  0^. 

FIGURE  8 

Two-color  TOFMS  of  toluene-propane  large  clusters.  The  energy  scale  is 

relative  to  the  Tol  0^.  The  concentration  of  propane  in  the  beam  is  2%  and  the 

nozzle  backing  conditions  are  P  =100  psia  and  T  =  4°C. 

0  0 


FIGURE  9 


Minimum  energy  configurations  of  Tol(C.^Hy)^  •  Atomic  positions  are 

given  in  Table  VII.  In  configuration  II  the  propane  interacts  less  with  the 

aromatic  cloud  than  in  configuration  I. 

FIGURE  10 

Minimum  energy  configurations  of  ToUC^Hq)^*  Only  configuration  IV  was 

y  o  d 

minimized  in  the  computer  modeling.  The  otner  configurations  are  composites 

of  the  geometries  of  Tol(C.,H„).  shown  in  Figure  9.  The  atomic  positions  of 

j  o  1 

configurations  IV  are  given  in  Table  VIII. 
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